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We analyze quark condensates and chiral (scalar) susceptibilities including isospin breaking effects 
at finite temperature T. These include rriu 7^ md contributions as well as electromagnetic (e 7^ 0) 
corrections, both treated in a consistent chiral lagrangian framework to leading order in SU{2) 
and 5'?7(3) Chiral Perturbation Theory, so that our predictions are model independent. The chiral 
restoration temperature extracted from (qq) = {uu + dd) is almost unaffected, while the isospin 
breaking order parameter {uu — dd) grows with T for the three-flavor case 5(7(3). We derive a 
sum rule relating the condensate ratio {qq) (e 7^ 0)/ {qq) (e = 0) with the scalar susceptibility differ- 
ence x{T) — x(0), directly measurable on the lattice. This sum rule is useful also for estimating 
condensate errors in staggered lattice analysis. Keeping m„ 7^ rrid allows to obtain the connected 
and disconnected contributions to the susceptibility, even in the isospin limit, whose temperature, 
mass and isospin breaking dependence we analyze in detail. The disconnected part grows linearly, 
diverging in the chiral (infrared) limit as T/M^, while the connected part shows a quadratic be- 
haviour, infrared regular as T^/Af^ and coming from vr^ry mixing terms. This smooth connected 
behaviour suggests that isospin breaking correlations are weaker than critical chiral ones near the 
transition temperature. We explore some consequences in connection with lattice data and their 
scaling properties, for which our present analysis for physical masses, i.e. beyond the chiral limit, 
provides a useful model-independent description for low and moderate temperatures. 
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The low-energy sector of QCD has been successfully described over recent years within the chiral lagrangian 
framework. Chiral Perturbation Theory (ChPT) is based on the spontaneously breaking of chiral symmetry 
SUL{Nf) X SUii{Nf) — SUv{Nf) with Nf = 2,3 light flavors and provides a consistent, systematic and model- 
independent scheme to calculate low-energy observables [TH3]. The effective ChPT lagrangian is constructed as an 
expansion of the form C = Cp2 + Cp4. + . . . where p denotes a meson energy scale compared to the chiral scale A^^, ~ 
1 GeV. For Nf = 3 case, the vector group symmetry is broken by the strange-light quark mass difference rris — rriu^d, 
although iTis can still be considered as a perturbation compared to A^^, leading to SU{3) ChPT, which reduces for- 
mally to SU{2) in the TOs — > 00 limit [3 . The formalism can also be extended to finite temperature T, in order to 
describe meson gases and their evolution towards chiral symmetry restoration for T below the critical temperature Tc 
[li 5^, where Tc ~ 180-200 MeV from lattice simulations [6,-^. The use of ChPT in this context is important in order 
to provide model-independent results for the evolution of the different observables with T, supporting the original 
predictions for chiral restoration [TU], also confirmed by lattice simulations, which are consistent with a crossover-like 
transition for Nf ~ 3 (2-1-1 fiavors in the physical case), which becomes of second order for Nf — 2, in the 0(4) 
universality class, and first order in the degenerate case of three equal flavors. 

The SUvi'2) vector group is the isospin symmetry, which is a very good approximation to Nature. However, 
there are several examples where isospin breaking corrections are phenomenologically relevant, such as sum rules for 
quark condensates [31 , meson masses jjll or pion scattering fT7, . For a recent review see [13] . The two possible 
sources of isospin breaking are the QCD — light quark mass difference and electromagnetic interactions. Both 
can be accommodated within the ChPT framework. The expected corrections from the first source are of order 
(m^ — m„)/m3 and are encoded in the quark mass matrix, generating also a tt'^t] mixing term in the SU{3) lagrangian 
[3]. The electromagnetic interactions are included in the ChPT effective lagrangian via the external source method 
and give rise to new terms pdUISllTS] of order L^^, Ce2p2 and so on, with e the electric charge. These terms are easily 
incorporated in the ChPT power counting scheme by considering formally — 0{p'^ /F'^), with F the pion decay 
constant in the chiral limit. 

The purpose of this paper is to study within ChPT isospin breaking effects related to the thermodynamics of 
the meson gas. We will be particularly interested in the physical quantities directly related to spontaneous chiral 
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symmetry breaking and its restoration, namely, the quark condensates and their corresponding susceptibihties at 
finite temperature. The quark condensate is the order parameter of chiral restoration, but since the transition is 
a smooth crossover for the physical case, different observables can yield different transition temperatures. Thus, 
the susceptibilities, defined as derivatives of the condensates with respect to the quark masses, provide also direct 
information about the transition and its nature, since they tend to peak around the transition point refiecting the 
growth of correlations. 

Let us mention some of the motivations we have in mind for the present analysis. For the physical values of quark 
and meson masses, we are interested in the effect of the isospin-breaking terms in the light quark condensate (uu + dd) 
and therefore on the ChPT estimates of the critical temperature. In addition, in the isospin asymmetric case, one 
has (uu) 7^ (dd) and in fact {uu — dd) can be considered an order parameter for isospin breaking. Actually, isospin 
is not spontaneously broken in QCD [16| which means that this order parameter should vanish for m„ = and 
e = 0. This is an important difference with the scalar condensate {uu + dd), which is nonzero in the chiral limit. It 
is relevant to estimate the thermal evolution of {uu — dd), since in principle the two condensates melt at different 
critical temperatures. A further motivation is the analysis of the three independent susceptibilities, directly related to 
the isosinglet, connected (isotriplet) and disconnected susceptibilities |T7j often discussed in lattice analysis [ T5H21) . 
Including properly the — vrid dependence of condensates is then essential to analyze the temperature and mass 
evolution of the connected and disconnected pieces measured in the lattice. In particular, the linear m^— m„ corrections 
to condensates survive the TOu = limit in the susceptibilities. The contributions coming from ■K^-q mixing in the 
SU{i) case belong to this type and are particularly important regarding the temperature dependence. This is not 
only interesting for physical masses but also to explore the scaling near the chiral limit, which in lattice studies has 
been used to investigate the nature of the transition |20| . In the lattice works, this scaling may be contaminated 
by lattice artifacts such as taste breaking in the staggered fermion formalism, which can generate contributions to 
susceptibilities masking the true scaling behaviour [201 121] . Our study provides then a model-independent setup for 
disentangling these effects and establishes the expected results in the continuum limit. 

We will work in ChPT to one loop, considering on the same footing the two sources of isospin-breaking. In a 
previous work ^TF\ we have studied the quark condensates at T = and several related phenomenological aspects of 
the isospin asymmetric case. We will refer to that work for more details about the formalism, the numerical values of 
the low-energy constants and other related issues. 

The paper is organized as follows. In section |ll] we will review the main aspects of the isospin-breaking ChPT 
formalism related to the present work. Our results for the quark condensates at finite T both in the SU{2) and SU{3) 
cases are given and analyzed in section [Till In that section we explore the temperature dependence of isospin breaking. 



as well as that of the sum rule relating condensate ratios. Section IV is devoted to the analysis of the different isospin- 
breaking scalar susceptibilities and their relation to the connected and disconnected ones. In subsection |IV A| we 
provide an interesting sum rule relating the electromagnetic differences in the condensates with the total susceptibility. 
We explore the possibility of using that sum rule to estimate the errors in the staggered fermion lattice analysis of 
the condensates, in connection with the taste breaking effect. In subsection IV B we make a thorough study of the 
connected and disconnected contributions to the susceptibility and their dependence with temperature, the quark 
mass and the isospin ratio mu/nid- We pay special attention to the connection of our results with different lattice 
analysis in the literature. 



II. FORMALISM 



The effective chiral lagrangian up to fourth order in p (a meson mass, momentum, temperature or derivative) 
including electromagnetic interactions proportional to is given schematically by Ceff — Cp2^^2 + £p4_|_g2j,2_(_e4 . The 
most general second order lagrangian is the familiar non-linear sigma model, including the gauge coupling of mesons 
to the electromagnetic field through the covariant derivative, plus an additional term proportional to a low-energy 
constant C compatible with the e 7^ symmetries of the QCD lagrangian pT] [T5|: 

Cp2+^2 = — tr [Df.U'^D^'U + 2BoM {U + U'')] + Ctr [QUQU^ . (1) 

Here, U{x) = cxp[i$/F] e SU{Nf), with $ the Goldstonc Boson (GB) matrix field for pions {Nj = 2) plus kaons 
and r] {Nf — 3), the latter being the octet member with I3 = S ^ 0. The covariant derivative is Z?^ = + iAfj_[Q, ■] 
with Af^ the EM field. Ai and Q are the quark mass and charge matrices, i.e., in SU{3) M = diag{mu,md,ms) and 
Q = (e/3)diag(2, —1, —1). Both the mass term and the charge one proportional to C in ([T]) break explicitly the chiral 
symmetry ^{/^(iV/) x SUuiNf) under which U — )■ LUK'^ with L,R G SU{Nf). The vector symmetry L = R is also 
broken for unequal quark masses and charges. Thus, in the light sector (w, d) the part of the mass term proportional 
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to m = (m„ + md)/'^^ the average light quark mass, is also proportional to the identity flavor matrix and therefore 
invariant under 5'C/y(2), while the part proportional to the mass difference ms = {niu — md)/^ and T3, the third 
isospin generator, is the one carrying out the QCD isospin breaking. The only remaining symmetry of the lagrangian 
^ is the U{1) L — R — exp{iXQ) corresponding to charge conservation. 

Working out the kinetic terms in ([T]) allows to relate the low-energy parameters F, BQmu^d,s, C to the leading-order 
tree level values for the decay constants and masses of the pseudo-Goldstone bosons. For SU{2) the masses read: 



In the SU{'i) case, the mass term in 
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induces a mixing between the 7r° and the 77 fields given by Cmix = 
{Bo/\/3){md — m„)7r''7y. This mixing between the two states with I3 = S = will play an important role in what 
follows. The kinetic term has then to be brought to the canonical form before identifying the GB masses, which can 
be easily done by the field rotation [3]: 



where the mixing angle is given by: 



7r° = TT^cose — T^sine, 
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Once the above n^rj rotation is performed, the SU{3) tree level meson masses to leading order read: 
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For pions, the main effect in the tt" — tt"'" mass difference comes from the EM contribution [23], while in the kaon and 

|24| (m„ = rud limit) indicate that m„ — md 



eta cases the violations of Dashen's theorem 



■Ml,=Ml^ 



corrections are relevant and must be kept on the same footing as the EM ones [Tl] . We emphasize that all the 
previous expressions hold for tree level LO masses with a = tt^, 7r°, K^,ri, in terms of which we will express all 
our results. They coincide with the physical masses to leading order in ChPT, i.e., M^^f^^^ = M^{1 + 0{M'^)) and 
so on for the meson decay constants = + ©(M^)). 

The fourth-order lagrangian consists of all possible terms compatible with the QCD symmetries to that order, 
including the EM ones. The C^i lagrangian is given in [5] for the SU{2) case, /ii,2,3 (contact terms) and Zi 7 denoting 
the dimensionless low-energy constants (LEG) multiplying each independent term, and in [3] for S'[/(3) the LEG 
named Hi 2 and Li ..iq. The electromagnetic Le'^p^ and Cf.i for SUl2) are given in O [13], A:i ...i3 denoting the 
corresponding LEG, and in 11 for 5'J7(3) with the Ki 17 LEG. The relevant terms needed for this work are given in 

m- 

The LEG are renormalized in such a way that they absorb all the one-loop ultraviolet divergences coming from 
Cp2 and £e2, according to the GhPT counting, rendering the observables finite and scale- independent. The numerical 
values of the LEG at a given scale can be fitted to meson experimental data, except the contact hi and Hi. The latter 
are needed for renormalization but cannot be directly measured, reflecting an ambiguity in the observables depending 
on them. The origin of this ambiguity is in the very same definition of the condensates in perturbation theory [5] . It 
is therefore convenient to define suitable combinations which are independent of those constants and therefore can be 
determined numerically. We will bear this in mind throughout this work and we will try to provide such combinations 
when isospin-breaking is included. The numerical values we will use for masses and low-energy constants in the S'[/(3) 
case are the same as in ^22] unless otherwise stated. In S't/(3) they come from the fits performed in [25] . 
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III. QUARK CONDENSATES AT FINITE TEMPERATURE 

The quark condensates for a given flavor qi at finite temperature T are given by: 



I d_ 



\ogZ 



dC 



drrii 



(6) 



where /3 = 1/T, V is the system volume, Z the partition function and {■)t denotes a thermal average. We will denote 
by {qq)T = [uu + dd)j' = —Jy]^ ^ogZ, the order parameter of chiral symmetry, while {uu — dd)T = ~Jv(j^ ^^gZ 
behaves as an order parameter of isospin breaking, since it is the expectation value of the part of the mass term in 
the QCD lagrangian proportional to to„ — to^ and e{qu — qd), respectively. It is still invariant under transformations 
in the third direction of isospin, which reflects electric charge conservation. 

In ChPT to one loop we obtain then the SU{2) finite temperature extension of the T = results in ||22 , which we 
give also for consistency: 



(99>'i 



{uu + dd)T - {qq)o + Bo [gi (M^o , T) + 2gi (Af,± ,T)]+0 (p^) 



{qq)o = -2F^Bo 



1 - - 2^^± 
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(7) 



where: 



(uu - dd)T = {uu - dd)Q = 'iBl{md - m„)/i3 - -P'^Boe'^kr + O [p^) , 

o 



(8) 



(9) 



and 



= 32^^°S7^' 
9^iM,T) . ^fdp^^-^^, (10) 

withE^ =p2_^M2. 

The expression Q contains the leading order tree level term from £2 given by {qq)o = —2F^Bo, the one-loop tadpole 
like contribution Gi{x = 0), with G the free meson thermal propagator, whose finite part yields the combinations 
fii+ gi{Mi,t) / (2F'^) (wc follow the same finite-T notation as in [5]) and the tree level from the fourth order lagrangian, 
showing up only at T = 0, which contains the LEG renormalized at the scale /z of dimensional regularization in the MS 
scheme [11[T3] so that the full expressions for the condensates are finite and scale-independent. Note that {qq)o includes 
the contact term hi. The isospin breaking in {qq)T for SU{2) is purely electromagnetic, showing up explicitly in the 

terms and implicitly through the pion mass differences. The temperature dependence is encoded in the functions 
5i(M, T) which increase with T and behave near the chiral limit (T > M) as gi{M,T) = ^^[1 + 0{M/T)]. 

Note that the effect of the electromagnetic corrections is to decrease the thermal part of {qq)T, since M^± > Mt^o. 
On the other hand, {qq)Q increases for the available estimates of the EM LEG, reflecting its ferromagnetic nature [22) . 
Our first conclusion is then that the critical temperature, estimated as that for which the condensate vanishes, increases 
with respect to the e = case, which is also a ferromagnetic-like behaviour induced by the explicit chiral symmetry 
breaking of the EM quark coupling in the QGD action. A simple estimate of the size of this effect can be obtained 
by taking the chiral limit m„ — uid = so that M^o = fi^o = 0, = 2Ce^/F^ and = VKPy^l + e^/Cj ~ 2/x^± , 
which gives T^'^^ /T^^^ ~ 1.003 with the parameters used in |22| and setting the involved ki to their maximum 
expected "natural" values k^ — l/(167r^). Thus, in principle we expect rather small corrections to chiral restoration 
from the electromagnetic breaking. Nevertheless, in section [IV A| we will go back to this point in connection with 
a sum rule relating the charge breaking with the susceptibility, suggesting larger corrections either for higher order 
transitions or for finite lattice spacing. 

The two sources of explicit isospin breaking in the lagrangian show up in the condensate difference ([s]), which 
depends linearly on to„ — uid with the contact /13 and vanishes for m,j = and e = in accordance with the 
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absence of spontaneous isospin breaking |TB| mentioned in the introduction. Recall that ^3 and fcy do not need to be 
renormalized and are therefore finite and scale-independent. An important point is that {uu — dd) does not receive 
pion loop corrections in the two-flavor case and it is therefore temperature independent to the one-loop order. In other 
words, isospin breaking in SU{2) does not change with T and the two condensates melt at the same temperature. 
This picture will change for Nf = 3 due to kaon loops and tt^t; mixing. 

In the SU{3) case, we calculate to one loop at finite temperature the light and strange condensates, taking into 
account both m„ — and e 7^ corrections. The condensates read now: 
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where: 



/c5_(m) 
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[6 (X7 + if8^(A^)) + 5 (X9-(m) + i^[o(M))] , 
[X9^(M)+i^[o(A^)], 



(14) 



In some of the above terms we have preferred to leave the results in terms of quark instead of meson masses. As in 
the SU{2) case, the results are finite and scale-independent, which concerns only the T — part [22]. 

There are some important differences with respect to the Nf = 2 case which deserve to be commented. First, the 
presence of the ir^r] mixing angle £ Q, as well as the more complicated dependence of meson masses with quark 
masses ([s]), imply that now m„ — corrections show up in {qq), apart from the EM ones. Note also that these 



corrections in {qq) and (ss) are at least C(£^) in the mixing angle, or equivalently in ' 



m^, except for an 0{e e) 



term in the kaon contribution. This is so because, apart from the explicit s dependence, one has to expand also the 
meson masses in ^ around £ = 0. All the masses depend quadratically on e except Af^± ~ —ae, Af^o ~ ae with 
a = {2Bo/ ^/i){mg — rh). Since, in addition, Af^± = Af^o + 2Ce^ / F"^ for £ = 0, we end up with the above mentioned 
term. 

Ano ther important difference between the two cases is that for S'C/(3) there are loop contributions to {uu — dd)T 
in (12). Kaon loops arise from the charged- neutral kaon mass difference, while neutral pion and eta ones from i^'^r] 
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mixing. When expanding in e now, t he l eading order is 0{s) even for e = 0. These hnear terms will be crucial for our 
analysis of susceptibilities in section IV Those loop corrections introduce now a T dependence in {uu — dd)T, unlike 
the SU{2) case. As it happened in the SU{2) case, we see that {uu — dd)T in (12) vanishes for and m.^ — nid, in 
agreement with 16], which we see from our analysis that holds including thermal corrections. 

At low and moderate temperatures gi{M.,^o,T) dominates over the kaon and eta contributions in (12 1, but it should 
be reminded that e in (|4| brings up a 1 /nis dependence which reduces the size of the pion term. In order to make a 
crude estimate, let us consider again the chiral limit, but keeping now the leading order in to„ — m^;, which we take 
then very small but nonzero while taking to — > 0+. In this limit the kaon masses are roughly kept to their physical 
values, which are well above the critical temperature. Thus, we consider the regime M.^- <C T <C Mk^ in which the 
pion term behaves as BQ{md — mu)T'^ /{^Anis) = BQ{md — mu)T^ / {18M^) . The kaon and eta contributions go like 

BqT^ [{rrid — to„)/tos] -y/ MK,ri/Te~^^'^'-'^^'^ [5], where we have taken also e = for simplicity. The pion term is still 
dominant due to the exponential suppression of K,r]. However, when compared to the T ~ part in that regime, 
which goes like {nid — TO„)i3g, we see that the quadratic growth with temperature is controlled by the scale 
instead of, say, the chiral restoring behaviour of {qq)T which is controlled by in the chiral limit. Therefore, the 
order parameter for isospin breaking {uu — ddjx grows with T, although it does so rather softly. Therefore, we do 
not expect big differences in the melting temperatures of the u and d condensates. This is also consistent with the 
expectation that in the limit where is arbitrarily large, say compared to to, the SU{2) result should be recovered, 
for which there is no temperature dependence for the condensate difference. 

The evolution with temperature of the condensate difference is shown in Figure [l] for the full case of finite pion 
mass and both e ^ and to„ 7^ rud- We have used the same set of low-energy constants and parameters as in [22], in 
particular rriu/md = 0.46 and ms/rh = 24. For the EM LEG Ki involved, we have displayed in the figure the two curves 
corresponding to their maximum and minimum expected natural values. We also show for comparison the result for 
TO„ — md, which shows that the charge contribution is actually of the same order as the one proportional to to„ — m^. 
We see that the T-dependent amplification of the isospin difference is rather large. In fact, this order parameter reaches 
values comparable to its T = value near the critical temperature, which is about Tc ~ 265 MeV in SU{3) ChPT. 
Nevertheless, due to the additional e suppressing factor discussed above, this enhancement is not enough to produce 
a sizable difference in the melting temperature of the u, d condensates, as it is clearly seen in Figure [l] (right), where 
we plot the two thermal condensates separately. The two plots showed in Figure [T] correspond then respectively to the 
two order parameters involved here: isospin breaking and chiral restoration. In turn, note that the curves on the right 
plot are independent of the choice of LEG since to this order {qiqiiT / {(iiqi)o = 1 — {{qiqi)T — {qiqi)o)/{BoF^) + 0{p'^) 
for i = u,d. 
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FIG. 1; Left: The u — d condensate difference (isospin breaking order parameter) at finite temperature in SU{?)), relative to 
its r = value. Right: The two condensates separately. 



The individual condensates in ( 11 1-( 13 ) contain the contact terms H2. These terms reflect an ambiguity in the quark 
condensates, inherent to their renormalization in QGD. It is therefore very important to deal with combinations of 
condensates which are free of this ambiguity. This very same source of ambiguity is also present in lattice simulations 
at finite T. A simple way to get rid of it is to subtract the T = contribution. This is the approach followed 
by the group [S] both for condensates and for susceptibilities. A different possibility is to consider the combination 
(99) ~ (jh) /■ms{ss) or for individual condensates in the isospin breaking case, {qiqi) — {mi/ms){ss) with i = u,d. 
Another sum rule free of contact ambiguities often used in T = phenomenology to relate condensate ratios [5| is the 
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following combination: 
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(15) 



where 0{mu — mdf , 0{e^) 0{e^{mu — rud)'^) have been neglected and Asr{0) is given in [52] with both sources of 
isospin breaking contribut ing at the same order, not only in the chiral counting but also numerically. 



We have seen in section 



III 



that the {dd)^ / {uu)t ratio receives significant corrections at finite temperature. On the 
other hand, we expect the strange condensate to vary slowly with T, from chiral symmetry breaking due to the strange 
quark mass. Therefore, we expect that the thermal corrections to this sum rule are also sizable. These corrections 



are plotted in Figure [tj for 
temperature. 



l/(87r ). They become comparable to the T = sum rule near the critical 
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FIG. 2: Thermal corrections to the sum rule relating condensate ratios. 



IV. SCALAR SUSCEPTIBILITIES AND ISOSPIN BREAKING 

In the isospin-breaking case, the scalar susceptibilities are defined as: 

Xij = -g:^{qjqj)T ^ 'py Qmdm " J J ^^^ii^i1^)i^''^)i^jqj)iO,0))T-l3V{q,qi)T{qjqj}T, hj^u,d,s 

(16) 

so that in the light sector, the relevant one concerning chiral restoration, we have three independent scalar suscepti- 
bilities Xuu, Xdd and Xud = Xdu- 

At this point, it is instructive to recall the definition of the connected and disconnected parts of the susceptibility. 
Consider the isospin limit with two light identical fiavors of mass m = m„ = rrid and e = 0. There is only one light 
susceptibility in this case, which can be written as: 



with: 
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XMs = ((Tri^r')')^ - (TrZ?r'>A, (18) 
Xcon - -{TtD-^)a. (19) 

where Di = — ih is the Dirac operator for every hght flavor in the QCD lagrangian and {■)a denotes integration 
over the gluon fields, so that formally Z — (exp^^ TrlogDj)^! where j runs over flavor and Tr runs over the space- 
time, Dirac and color indices. This separation is important for lattice analysis, as we will discuss below, and reflects 
the contributions with connected and disconnected quark lines, since D^^ is the quark propagator. However, when 
considering the low-energy representation for the partition function and the susceptibilities in terms of GB fields, it is 
not so simple to separate the connected and disconnected parts if we take the isospin limit from the very beginning. A 
possible approach to perform such separation is to work within the partially quenched ChPT framework, as discussed 
in [27j for the vacuum polarization. We will however work within the isospin breaking scenario we are considering 
here, which is very useful for this purpose, as noted first in [T7] for the susceptibilities and used also in [55] for the 
vacuum polarization. The main point is that for ^ md' 

Xud = ( (Tri?- 1 ) {TrD- ^ ) ) ^ - {TrD^ ' ) a (Tri?- a, (20) 



so that one has Xdis = linim„^md X«d and since dm — 9m„ -I- dm^, from (17) we ha^ve cilso \con — 

limm^^ma [iXuu + Xdd)/^ - Xud]- 

Therefore, with this observation in mind, we define in the isospin-breaking regime the following basis of total. 



connected and disconnected susceptibilities in terms of the ij basis in (16): 



X = Xuu + Xdd + ^Xud, (21) 

Xcon — 2 iXuu Xdd) Xudj (22) 
Xdis = Xud, (23) 

which we can therefore obtain directly from our expressions for the isospin-breaking condensates obtained in the 
previous section. Observe that none of the Ki dependent terms in the condensates depends on the quark masses and 
therefore the susceptibilities are independent of the EM LEG. 



Note that according to (16), x in (21) corresponds to the correlator of the isosinglet condensate (qq), the order 
parameter of chiral restoration, while the connected contribution Xcon is the correlator of the isotriplet uu — dd, the 
order parameter for isospin symmetry. A divergence or sudden growth of these susceptibilities would indicate then a 
phase transition for the corresponding order parameter. 

We also remark that the definitions of the connected and disconnected parts in terms of uu, dd, ud ones are not 
unique. We could as well have defined Xdis as a(Xim~Xdd) + Xtid for arbitrary a, which also reduces to the combination 
(18) in the isospin limit. We are following the same convention as jl7]. These formulas can be easily extended to Nj 



identical flavors, for which x = NfXcon + Njxdi 

In the following we will analyze several aspects related to the above defined susceptibilities in different limits. 



A. Sum rule for EM-like corrections to condensates 



Before studying in detail the different susceptibilities, in this subsection we will relate the EM corrections (and 



actually any charge-like correction to pion masses) to the condensates, found in section III with the total scalar 
susceptibility. Gonsider first the condensate calculated in SU{2) in ([t]) and let us define the ratio: 

Now note that to one loop, the explicit dependence of the condensate in is only in the T = part, since the 
charge dependence in £2 is contained implicitly in the pion mass differences. Therefore, r{T) — r(0) depends on the 
charge only through the parameter 5^^ = (Af^± — M^o) /M^o, in which we can further expand (for the EM pion mass 
difference S^^ ~ 0.1). Taking also into account that the condensate is just the sum of the tadpole contributions for the 
three pions, we can write: 
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r(T) - r(0) 



Ml 



d 



Ml d_ 
6B'^F^ ""dm 



[{qq)T~{qq)o]+0{6l) + 0{p'), 



which, from the susceptibihty definition in (21 1 can be written, to this order, as 



(25) 
(26) 



r{T) - r(0) = 



3 M2F2 



[X{T) - X(0)] 



(27) 



This sum rule relates then pion mass deviations in the condensate with the total scalar susceptibility. Note that the 
above result is written only in terms of the quark mass, the pion mass and decay constant and the charged-neutral 
mass difference, without specifying if the latter is of electromagnetic origin. It states that, even though the mass 
deviation J.^- may be small, the corrections to the condensate may be amplified near the phase transition, where the 
susceptibility is maximum, if such transition is sufficiently strong. Actually, the quantity proportional to St^ on the 
right hand side of (271 is directly measurable on the lattice [8^. 

For the case for the electromagnetic mass difference in SU{2) discussed in section |llj we have i5^A/^o = 2Ce^ jF"^ 
and: 



r{T) - r(0) = — g2(M,o, T) + ©(e^). 



(28) 



with: 



r(0) 



4Ce2 
F4 



(29) 



and 



F 



1 



92{M,T) 



d9i{M,T) 
dM^ 



327r2 

_ 1 



1 + log 



dp 



/i2 

1 



1 



Er, ef^^p - 1 ■ 



(30) 



Note that r(T) is finite, scale- independent and also independent of the e = LEG l^,hi,h^. In particular, it is free 
of the contact-terms ambiguity, which makes it a quantity suitable for physical predictions. It is also independent of 
-Bo, unlike the individual quark condensates, which have only physical meaning when multiplied by the appropriate 
quark masses, since the miBo products give meson masses. In addition, the dependence with the EM LEG disappears 
in the difference r{T) — r(0), which is the quantity directly related to the susceptibility through (27). 

The above relation can also be explored for S'[/(3). However, the connection with the susceptibility is not dire c t in 
that case. The dependence of {qq)T— {qq)o enters now through M^± and M/f ± . However, the condensates ( 11 1-( 13 ) 
depend on the light quark mass through all the meson masses M^,Mjf,M^. The result is that r{T) — r(0)can be 
expressed as the susceptibility term in (27) plus a linear combination of d{{ss)T— {ss)o)/dm and 9((ss)t — (ss)o)/9m 



to this chiral order and neglecting 0{e'^) and Oiniu — irid) isospin-breaking corrections in the right hand side. Since 
the strange quark condensate has a much weaker dependence on temperature than the light one (or equivalently, we 
can approximately neglect the thermal functions evaluated on kaon and eta masses) we expect the T behaviour of 
r(T) — r(0) to be dominated by the light scalar susceptibility also in the SU{3) case and therefore the sum rule (27) 
should hold approximately. In this case we have to one loop: 



r(T)^^(3)_^(0)Sa(3) ^ l + ^[2g2(M,±,r)+.g2(M^±,T)]+0(e4), 



r(0) 



5(7(3) 



l + e2/CS+(M) 



2Ce2 



where the expansion in to leading order allows to express the result in terms of the tt^ and masses. As in 
SU{2), r{T) is finite, scale-independent and independent of the e = LEG, so that it is free of contact ambiguities. 
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We compare the above expression with the susceptibihty in the SU{3) case. As it will become clear in section IV B 



the corrections to the total susceptibility x from both sources of isospin breaking are small. Actually, we will see that 
the NLO correction in the QCD breaking is ©(m^ — mj)'^. The one-loop result is: 



X(T) - x(0) 4m2 



^0 



M4 

x(0) 4™^ 



352(M., T) + g2{MK,T) + -g2(M^, T) 



[X(T) - x(0)] = 2 
X(0) = 16 [8Ll{p) + 2Ll{^J) + - 4 



2,V^ + VK+ - Vq 



(31) 

2 



(32) 



The two quantities are compared in Figure [3j The deviations between them are negligible for the range of rel- 
evant temperatures. Therefore, although for physical masses the electromagnetic corrections are relatively small, 
they grow with the susceptibility, which is a model-independent prediction. For comparison, taking the value 
of {rh?) / {M^) [x{Tc) — x{^)] from the lattice simulations in [S] for 2-fl flavors with the lattice Tc value gives 
r(Tc) — r(0) ~ 0.013, not far from the higher temperature values in Figure |3j although the ChPT curve cannot 
reproduce the susceptibility peak, only the low and moderate T behaviour. These small EM corrections for the con- 
densate are in accordance with our simple estimates made in section |III| and translate into a few MeV difference in 
the determination of the critical temperature from the order parameter. 




150 
T (MeV) 

FIG. 3: The r{T) function in SU{3) encoding the EM corrections of the quark condensate. We compare it with the normalized 
light susceptibility for physical quark and meson masses. For these values, r(0) ~ 1.01. 



The sum rule (27) has another interesting consequence, regarding lattice simulations. In the staggered fermion 



lattice formalism, the need to introduce four different copies (tastes) for every quark flavor leads to the so called taste 
violation [HI [191 ES] . This is a lattice artifact which in some aspects is similar to the isospin or flavor violations we 
are analyzing here. The new tastes enlarge the chiral symmetry group to SU{4Nf) x SU{4:Nf), producing then 15 
pseudo-Goldstone bosons plus one massive state (r^'-like) for every quark flavor. All these new meson states become 
degenerate in the continuum limit, where taking the fourth root of the Dirac fermion determinant is enough to remove 
all the spurious copies. However, for finite lattice spacing a, the tree-level masses of those states receive O(a^) 
contributions, which break explicitly the chiral group in the lagrangian, only one Goldstone boson remaining massless 
in the chiral limit leaving then a residual 0(2) or C/(l) symmetry. The mechanism is similar to the electric charge one 
we are analyzing here, by which the charged states receive 0{e^) corrections and the U{1) EM symmetry remains. In 
fact, for the staggered case one can construct a generalized chiral lagrangian including all possible terms compatible 
with the new symmetry. This is called staggered Chiral Perturbation Theory |29l 130] . Among the new terms one 
recognizes contributions of the form Tr[^C/^C/^] with ^ a given combination of SU{A) generators, i.e., like the charge 
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term in ([I]) in 5*17(3). Obviously, the staggered case includes other operator structure and the spectrum of states 
is more complicated. However, we can use the sum rule (27 1 to estimate roughly the expected differences between 
the lattice staggered condensate and the continuum one, considering the lightest states. For fine enough lattices, 
one has 5° = (M^^^ — M^)/AI^ ~ ca^ for the lightest tastes of squared mass M^^^ [5], where from the two smallest 
lattices in [B] we get c ~ 140 fm^^. With this S'^ we can then use (27) to estimate r°'{T) — r°(0). Consequently, we 
expect the larger errors coming from this taste violation effect to appear near Tc- That is indeed the case when we 
compare lattices of decreasing temporal extent Nt — a/T for the condensate data given in [8 . More quantitatively, 
taking also the susceptibility values of [S], we get r°'{Tc) — r°(0) ~ 0.07ca^. Estimating the T = part using (29) 
with — > S'^M^F'^ / {2C), we get a relative correction for the condensates near Tc with respect to the continuum of 
about 20% for the Nt = 12 data in [8] and about 12% for the Nt — 16 ones in [3T]. Following the same idea, we 
get a relative difference between the Nt = 10 and Nt — 12 lattices of around 8% near T^, which is actually in good 
agreement with lattice data [S|. A direct translation into an error for the critical temperature is not easy to obta in. 
If we simply extrapolate the one-loop chiral limit expression {qq)T — (99)0(1 — T'^/T^), writing the l.h.s of (27) in 
the one-loop equivalent form (25), we get a very rough estimate ATc ~ 10 MeV for Nt = 12 with respect to the 
continuum, although the chiral limit is not always numerically accurate, as we will actually see in the next section. 

Estimating taste violation effects is important, since they are one of the main sources of the discrepancies between 
different lattice groups for the determination of the critical temperature. An important effort has been made over 
recent years to minimize these effects, not only by considering finer lattices, but also by introducing lattice actions 
where taste symmetry is reduced [5T] . 



B. Temperature and mass dependence of connected and disconnected susceptibilities. Relation with chiral 

restoration and lattice analysis 

The behaviour of condensates and susceptibilities with temperature and quark masses is crucial in order to un- 
derstand the nature of the chiral phase transition when approaching the chiral region (mg,T) — )■ (0+,Tc). Lattice 
simulations have addressed the question of how those quantities scale with and T until very recently |20j . An 
essential part of this program concerns the scaling of the connected and disconnected parts of the scalar susceptibility. 
The disconnected piece is given in terms of closed quark lines and is therefore directly related to {qq) and expected to 
be sensitive to chiral restoration. Actually, near the chiral limit, i.e., the infrared (IR) behaviour, it is known to scale 
as Xdis ^ ^'^'w for T = and Xdfs ^ T/Mt^ at finite temperature 117 . The infrared contribution is then controlled 
by the GB loop contributions. The situation is not so clear for the connected part, since its infrared divergent piece is 
proportional to nj: — 4, with Uf the number of identical light flavors |17j and therefore vanishes for n/ = 2. However, 
its IR finite part contributes in the physical case of massive pious and is actually an important difference between 
QCD and 0{N) models in which the lattice scaling fits are based Besides, the connected contribution receives 
important "false" GB-like corrections coming from taste violation 20, 21]. It is therefore important for lattice studies 
to provide the continuum result for the disconnected and connected susceptibilities in the physical case of 2+1 flavors 
and massive pions. 

Our present ChPT one-loop analysis allows to obtain a model-independent prediction for the low temperature and 
small mass behaviour of the susceptibilities. The inclusion of isospin breaking effects is crucial. In fact, it will be 
useful for the following discussion to note that: 

Xco7i = ^ iXuu - Xud) + ^ {Xdd - Xud) = {uu - dd), 

Xdis = -7 [drh{uu + dd) - 2dms - dd)] , (33) 



with mg ~ {rriu — md)/2, so that Xcon comes only from the condensate difference in (12) and its leading order is 
obtained from the linear terms in m„ — m^. 

Now, from our considerations in section |IIIt if we neglect for the moment the charge corrections and we take the 



quark mass derivatives in (11) we have, to leading order in ms, Xuu — Xdd, Xuu + Xud — x/2 and Xuu - Xud — X 



con 



where the leading terms Xi Xcon and Xdis = x/4 ~ Xconl'2 are 0(1) in the m,? counting and are given in S'C/(3) by 



(31|-(32| and 
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1 



3(M2 - Af2) 



Xcon Xcon (0) 



R2 
^0 



92{Mk.T) + 



2 [gi (M^, T) -.91 (M^,T)] 2n , ( md - rriu ; , 2x 



(0) 



3(M2 - M2) 

(^^ - ^,,) 



R2 



(34) 
(35) 

(36) 

(37) 



Several remarks are in order from the previous expressions. First, we emphasize that all of them are finite and 
scale-independent, which can be explicitly checked from the scale dependence of the LEG [Hj. We also note that, 
unlike the disconnected part, the connected susceptibility does not receive contributions from the mass derivative of 
pion tadpoles, i.e or 52(-^7ri T). These turn out to be the dominant ones in the chiral limit (see below) and this is 
what we expected from our previous d iscu ssion on the infrared behaviour. Thus, we identify the difference /x^ — /i^ in 
(37) and the corresponding gi one in (36) as the contribution of tt^t] mixing, while the vk , g2{MK) terms come from 
the expansion of the kaon contribution in the r.h.s. of ( 12 ) around the isospin limit. From ( |33[ ), the disconnected part 
receives in addition a contribution from the sum (11 ) and hence it incorporates the z^^, g2{M^ 
the T — part of Xcon depends on the contact term H2, while that dependence cancels in Xdis 
quantities such as Xcon{T) — Xcon(O) can be unambiguously determined, similarly as the quark condensate case (see 
our discussion at the end of section III). This is a relevant comment for lattice evaluations of this quantity. In fact, 
as a consequence of the vanishing pion terms, we see that Xcon{T) — Xcon(O) vanishes formally in the mg — > 00 limit, 
recovering the pure SU{2) result that we would get from ([s]), which holds also for the other susceptibilities taking 
into account the conversion between the SU{2) and SU{3) LEG [52] (see below). 

Regarding the expansion in (m^ — m^) , it can be seen from (|33[) , the condensate expressions (111, ( 12 ) and the meson 



terms. Note also that 
This means that only 



masses and mixing angle dependence on rh and ms, that the linear order in mg cancels both in the connected and the 
disconnected parts and so we have written in the previous expressions and in the total susceptibility (31 1- (32). It is 
important to remark that the 0{ms/ins)'^ ^ C'l^)^ corrections contain "tadpole mass derivative" terms ^^,(77(M^) in 
both 

Xcon s-rid Xdis- Howcvcr, an important difference between them is that those IR dominant terms do not appear 
to leading order in the connected contribution. In fact, since v^o = -\- [dv-^ j dM'^)0{ms / , where the subscript 
TT indicates just the m„ = rud and — Q pion, the disconnected part receives an 0{ms /msY' proportional to the 
second derivative of the pion tadpole and hence more IR divergent, and so on for the thermal part g2{Mj^o, T). This 
contribution is not present in the connected part to that order. Thus, we expect the isospin-breaking corrections to 
be larger in the disconnected than in the connected susceptibility. 

Another pertinent comment is that we have been able to obtain the leading order for Xcon and Xdis only after 
considering properly all the m„ 7^ md contributions and then taking the m„ = nid limit. However, one can be led 
to misleading results by setting to„ = from the very beginning. For instance, for two equal masses one could 
think naively that Xuu — Xud — Xdd = x/4, from the definition (21). However, from our previous analysis we see 
that in the isospin limit what we get actually is Xud = x/4 - Xcon /2and Xuu = Xdd = x/4 + Xcon 1 2 with Xcon given 
in (36)-(37) i.e, not vanishing for m„ = and physical m^, although formally suppressed in the 00 limit. 



Thus, in the isospin limit taking just one flavor susceptibility and multiplying by four does not give the total scalar 
susceptibility^ which should be obtained instead by considering the derivative of the full sum of u and d condensates 
as given by (21). Note that this correction does not affect to condensates, for which the {uu) and {dd) difference 
vanishes for niu = rud and = and therefore {qq) = 2{uu) — 2{dd) in the isospin limit. This comment may be 
relevant for certain lattice analysis, where working within the "one-flavor" equivalent framework is often done, since 
in that way it is easy to discuss for instance the taste-breaking effect. This observation could then help to explain 
the worse 0{N) scaling properties of the susceptibility with respect to those of the condensate [20]. In that work, 
the lattice data for the susceptibility scaling function in the 2-1-1 case suffer from a sizable increase as the strange 
quark mass is decreased relative to the light one. That increase is not seen in the quark condensate data and could 
be due partly to the deflnition used as we have just explained, since the positive term proportional to Xcon increases 
with l/nis, from ( |36| -(37). This is not the only effect that may cause this "wrong" scaling, since, as pointed out in 
[201 121| taste-breaking induces an artificial infrared pion contribution in Xcon which is not present in the continuum, 
as our above expressions show. What we are pointing out here is that the "one-flavor" Xuu and Xdd are sensitive to 
isospin-breaking terms even for = nid and in the continuum, unlike considering for instance the total susceptibility 
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X ^ Xuu 

by its T 



'Xud in (21 1 for which the Xcon term cancels. In addition, as we will see below, Xcon is dominated numerically 



part, which would explain why the anomalous scaling is reduced for the subtracted susceptibility. 
The charge corrections to susceptibilities in our previous expressions ([3T|)-([37| arise only through the mass of the 
charged mesons and . 
and 2iyK — >■ vk - 
and 2g2{MK,T) 



ace 



52 (M, 



in ([32l) and ^ 



Thus, including the charge amounts to rep ; 
in ( [32| ) and (37) and so on for the thermal parts 3g2{MT^,T) — > g2{M^a,T) + 2g2{M^±'^T) 
-> g2(Afic,T) + g2{Mji± ,T). Although for physical values of the electric charge and masses, these 
represent small perturbative corrections, the fact that near the chiral limit the coefficient of the IR-dominant and 
92{MTr,T) reduces in 1/3 for 7^ is the reflection in the scalar susceptibility of the behaviour of the condensate in 
terms of corrections to the masses analyzed in section IV A Thus, when the mass corrections (5^ become sizable, 
as in the staggered lattice formalism, the susceptibility is reduced by that factor, which eventually would imply that 
the transition peak or maximum is displaced to a higher temperature, consistent with our analysis in the previous 
section about the increasing of in the condensates. Recall that one cannot just expand the pion terms in and 
then take the chiral limit, since that expansion assumes that the charge part of M^± is small compared to the quark 
mass one. 

Before continuing, we also remark that our above results are compatible with the recent observation [2T that the 
connected and disconnected susceptibilities can be inferred from the zero momentum limit of the ap and /o correlators 
calculated previously in staggered ChPT [32] . The motivation of those works is precisely to estimate the contribution 
of heavy pion-like tastes to the IR part of Xcom which could mask the scaling behaviour. The continuum limit of the 
results in plj reveals the same tt, if, 77 loop contributions as in our expressions (34|-(37). We provide the full ChPT 
result, including the LEG contribution necessary to guarantee the finiteness and scale-independency of the results, as 
well as the analysis of the higher order corrections in isospin breaking. 

As discussed above, the behaviour of the susceptibilities near the chiral limit (IR regime) is very illuminating 
regarding their approach to chiral restoration within the 0(4) or SU{2) pattern in the continuum, i.e., without taste 
breaking effe cts. Let us consider this regime first for T = and only for the leading order terms in the isospin 
expansion in (35l-(36l, i.e., we set to,j = = m and = in the T — Q susceptibilities and consider m ^ rris. We 
denote by a superscript "IR" the nonvanishing terms in that limit: 





"^"'J. °' = -5!;^ '"S 7? + 32i6(rt + I -28 + 51og I , (38) 



The IR divergent log term in (38 1 coincides with the one obtained in [T7j, where a cutoff regularization was used. 
Multiplied by 4, this is also the IR dominant part in the total susceptibility x- In addition to that term, we obtain here 
the regular part, not IR divergent but not vanishing in the chiral limit, which provides the dependence with the LEG 
and together with the logarithm gives the consistent scale-independent GhPT prediction. As for the connected part, 
again only the IR divergent contribution is given in |17j , which as commented above turns out to vanish exactly for two 
light fcrmions. Here we also give the regular contribution, also scale-independent, which unlike the disconnected part 
depends on the contact term iJj- As an interesting consistency check, we can recover the SU(2) limit from the above 
expressions, using the conversion between the LEG of SU{2) and SU{?>) [3 + — — 16Lg + 5i^,,/18 — l/(967r^) 

and hi = ALl + 2H^ - i^k/^ - v-^l'^ + l/(967r2), so that x^^^^^ = -"iv^ + 2{ll + h\- h^) and Xcon^^ = 4/1^ , which is 
the same result that we would have starting directly from the SU (2) expressions in ^ and ^. 

Let us consider now the dominant IR thermal contribution in the 2+1 flavor case, i.e, apart from rh <C mg we 
also consider temperatures ^ T ^ Mk, so that we neglect all the Boltzmann exponentials exp(— Mif^^/T) and 
expand gi{M^,T) = f^[l - 3M^/(7rr) ©(M^ logM^)] and g2{M^,T) = T / {%t:M^) + 0{\ogMl) [51. Thus, we get: 

[XdtsjT) - Xdj^(O)]^^ _ 3r 

Bl 167rM/ ^ ' 

\Xcon(T)-Xcon{Q)r _ 



(41) 



The disconnected part (40 1 is again the one obtained in [T7| in the IR limit. It diverges more strongly than the 



T = contribution in (38) in this limit, revealing its critical behaviour. The growth with T is linear over the GB 



mass scale. Recall that, apart from the heavy masses thermal exponentials, we are neglecting also logM^ terms in the 
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disconnected part ( |40[ ). The situation is completely different for the connected contribution, which is regular, albeit 
not vanishing, in the chiral limit. The quadratically growing term in (41 1 survives for — >■ against neglected 
©(M^) and is dominant over exp{—Mfc,^/T). It vanishes formally as — > oo, recovering the SU{2) limit. For 
physical masses though, an specific and model-independent difference between the Nf — 2 and Nf = 2 + 1 cases is 
the (soft) temperature dependence of the connected susceptibility, the scale that controls its growth being instead 
of the of the connected part. This is a consequence of Xdis measuring the fluctuations of the chiral restoration 
order p ara meter, while Xcon is related to those of the isospin-breaking one, i.e., {uu — dd), which as we have seen in 
section III increases moderately. Note that near the chiral limit we could as well have written the term divided by 
just by changing the multiplying factor. Keeping reminds its vr'^?/ origin, as it is clearly seen in the original 
expressions in ([36|-(37). 



In Figure[4]we plot our numerical ChPT results for the susceptibilities, including all the isospin-breaking corrections. 
The plots in that figure show the difference with respect to the T = results, which are collected in Table |Ij At T = 
we use the same LEG values as in [52], which are quoted in the table. Remember that the susceptibilities are 
independent of the EM LEG and that the disconnected one is independent of contact terms. The contact LEG H2 
appearing in the connected contribution is estimated from resonance saturation arguments |151l26j . The normalization 
used Bq = Af^/(4?fi^) is the same one used in some lattice works [S]. 





Xdis /Bo 


Xcoti/Bq 


x/B'i 


4xun/B'i 


ms/rh — 24 


0.024 


0.025 


0.146 


0.196 


nis/m = 24, IL 


0.025 


0.025 


0.148 


0.197 


m,s/rh = 10 


0.016 


0.024 


0.113 


0.163 


ms/rh = 10, IL 


0.017 


0.023 


0.114 


0.161 


ms/rh = 100 


0.036 


0.026 


0.194 


0.245 


ms/m — 100, IL 


0.038 


0.025 


0.203 


0.254 



TABLE I: T = values for the different susceptibilities in ChPT. The isospin limit (IL) values correspond to m„ — > md and 
e — 0. For the third to sixth rows we fix m^/md and and vary the light to heavy quark mass ratio. The first and second 
rows correspond to the physical values. The LEG values used are = 2LS = 1.24 X 10-^ Lj; = at the scale fi = 770 MeV. 



We see in the plots that the general features explored in our previous ana lytical discussion are well reproduced. 
First, the C(e^) and ©(m^j — to„)^ terms neglected in (31)-(32| and (34)-(37) are numerically small for the relevant 



temperature range, for physical values of quark and meson masses. In fact, as anticipated in our previous discussion, 
we see that those isospin corrections are larger for the disconnected than the for connected part and are also larger 
with temperature, all due to the appearance of IR terms proportional to the second derivative of the tadpole in 
Xdis- Remember that the leading order in the isospin limit comes actually from the 0{mj^ — niu) terms in the 
condensates. Second, we appreciate qualitatively the linear and quadratic growth with temperature of the disconnected 
and connected parts respectively, as expected from the infrared analysis. In fact, we see that although the connected 
term grows faster, its absolute value is much smaller due to the scale compared to the of the disconnected 
part. However, it is important to note that the IR limit expressions (40)-(41) are numerically rather far from the 
exact ones for the physical pion mass. The difference is larger for the disconnected contribution, since, as stated 
above, in (40) we are neglecting ©(logM^r) terms, while in (41 ) the neglected terms are O(Af^). In fact, this justifies 



further our present analysis, since we provide the full expressions beyond the chiral limit. Also as discussed above, 
the infrared limit expressions for T = given in (38l-(39) survive not only the chiral limit but also the — >■ 00 
one, Xdis still diverging but only logarithmically, which for physical masses makes the two susceptibilities numerically 
comparable. This is clearly seen in the values given in the first two rows of Table |l] Actually, for this very same 
reason, and following our previous discussion, the deviations of Xuu from the naive isospin-limit expectation x/4 is 
much more pronounced at T = than for finite T, as it can be seen by comparing the last two columns in Table [T] 
which give about a 30% relative difference, while the last plot in Figure [4] where we compare their thermal differences 
give only corrections below 10%. In fact, this is consistent with our previous discussion about the influence of the 
connected part in the scaling properties observed in the lattice. If we consider the subtracted susceptibilities as defined 
in |20J from the subtracted condensate {uu) — (mu/mg){ss) , we see that the dependence on the LEG disappears in 
the rh/rria — > limit. Remember that in this limit all the T = contribution of the connected part is absorbed in 
/13 (see our previous comments) and therefore considering the subtracted susceptibility is equivalent to switch off the 
dominant T = part of the connected susceptibility. This is indeed observed in the lattice [5D] since the subtracted 
susceptibility fits better the expected 0{N) scaling behaviour than the unsubtracted one. 

The variation with the quark mass is displayed for T ^ in Figure [5] and for T = in Table |lj It is important to 
remark that we have chosen to keep fixed the ratio niu/md — 0.46 (same value used in our previous analysis [22j ) and 
m<,, while we vary m/mg above and below the physical quark mass ratio. In other words, when m — > 0"^, — > 0"*" 
while Mf^ri remain fixed. This is meant to be the relevant limit when approaching chiral restoration. In addition. 
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FIG. 4: Temperature dependence of disconnected, connected and total susceptibilities in three-flavor ChPT, for physical quark 
and meson masses. The full blue curves show the results with all the isospin-breaking corrections of higher order included, while 
in the dashed red ones we display the leading order in the isospin limit. We also show the IR expressions for the connected and 
disconnected parts, as well as the deviations of 4xuu from xtot in the isospin limit (last plot). The T = Q results are given in 
Table m 



smce we can write {md — mu)/ms — 2{rh/ms){l — mu/md){l + mu/md)^^ , s scales in this limit as 0{rh). Recall 
that, although the values of the are fitted to low-energy data with physical masses those LEG are formally 
independent of the quark masses [3j. The same applies to the tree level value of F we are using. 

As we expected from our previous IR analysis in eqns. (40) and (41 1, the light quark mass dependence of the 
thermal disconnected susceptibility is much stronger than the connected one, as seen clearly in Figj5] and as long as 
we take the limit in the order specified above. In terms of chiral restoration, this anticipates a much stronger growth 
or peak near Tc for the disconnected part. From the same arguments, the behaviour of the connected part is expected 
to be softer near the transition, although growing with for low and moderate temperatures. We also show in the 
figure the comparison with the infrared limit for the smaller rh/rrig case, where it can be seen that the curves are now 
closer than for the physical pion mass case in Figure |4] 

In addition, the isospin corrections are also more important for the disconnected part, where they actually increase as 
TO is decreased, than for the connected one, where the isospin limit and complete curves are almost indistinguishable 
in the figure. The same holds for the T — contributions in Table [1] According to our previous discussion, this 
behaviour of the isospin corrections arises from the dominant IR terms Xdis ^ B'^{T M"^ / M^)e'^ = 0{^/m) as compared 
to Xcon ~ Bq{T / M-,^)e'^ — ©(to'^/^). This effect is weaker for the T — contributions since the IR leading corrections 
i^TT^^ (in Xcon) and {dv^^ / dM'^)e'^ (in Xdis) diverge softly, namely, as m^logiTL and to respectively. Thus, although the 
isospin corrections are amplified in the disconnected susceptibility for large temperatures and small masses, they are 
still perturbatively under control in the chiral limit. 

The limit where to/tos vanishes not by taking to — >■ 0"*" but keeping to fixed and taking TOj — cxi is where we 
recover the pure SU{2) results, as discussed before. In that case, it is the connected part which is more sensitive 
to the quark mass variation, vanishing for large Af^, while the disconnected one remains invariant. Although this is 
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FIG. 5: Quark mass dependence of the thermal disconnected, connected and total susceptibilities in three-flavor ChPT, for 
fixed TUs and fixed rriu/md- We also show the isospin limit (IL) m„ — >■ md and = in the smaller light mass case. The 
T = results are given in Table |l] For comparison, we also display the infrared limit in the ms/rh = 100 case. 



formally interesting for connecting the SU{2) and SU{3) cases, it is not so relevant for studying the critical behaviour. 



V. CONCLUSIONS 



In this work we have analyzed the relevant observables regarding chiral symmetry restoration, namely quark con- 
densates and scalar susceptibilities, in the presence of isospin breaking. We have considered on the same footing the 
QCD {mjj^/m^ mass difference) and electromagnetic corrections, to one loop in Chiral Perturbation Theory, both in 
the SU{2) and SU{3) sectors. Our analysis provides useful and model-independent results regarding several relevant 
aspects of isospin breaking and chiral restoration, which may be particularly interesting for lattice studies. 

The sum {uu + dd)T, the order parameter for chiral restoration, receives small isospin-breaking corrections for the 
physical values of masses and electric charge. These corrections affect only slightly, less than 1%, to the value of the 
critical temperature, which they increase as a ferromagnetic response. The difference {uu — dd)T is the order parameter 
of isospin breaking. It is temperature independent in the SU(2) limit, but when kaons and eta are included, it shows 
an increasing behaviour, which in the chiral limit is given by (m„ — rriiijT'^ / . The deviations with respect to its 
T — value become sizable as the temperature is increased, but they are controlled by a larger energy scale 
than the typical of {uu + dd)T- This large growth of isospin breaking does not reflect in the chiral restoration 
temperatures of {uu)t and {dd)T, which remain close to each other, consistently with the idea that chiral restoration 
is little affected. We have also evaluated the temperature corrections to the sum rule relating the {ss)t/ {uu)x and 
{uu)t I {dd)T ratios, which is useful because it does not involve undetermined contact low-energy constants. The 
corrections in this case come directly from the {ddjx/ {uu)t ratio and are therefore rather large for the temperatures 
of interest. 
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A very important part of the present work has been the analysis of scalar susceptibilities in the isospin asymmetric 
scenario. We have related the different flavor susceptibilities with the total, quark connected and quark disconnected 
susceptibilities often used in lattice analysis. Electromagnetic corrections to the quark condensate turn out to be 
directly related by a sum rule to the total susceptibility and then to the growth of fluctuations, which is meant to be 
maximum near the critical point. This sum rule is valid for any small deviation of the pion masses, as for instance the 
one arising in the staggered lattice formalism due to taste breaking effects. Actually, we have made rough estimates 
of the corrections to condensates expected from this source, comparing lattices of different sizes among them and with 
the continuum limit. These estimates are in good agreement with the errors quoted in the lattice works. 

The isospin asymmetric calculation allows for a direct extraction of the connected and disconnected susceptibilities, 
even in the isospin symmetric limit. The terms in {uu — dd)T linearly proportional to rriu — rnd give contributions to the 
connected part not vanishing in the isospin limit and which affect for instance the naive extrapolation of a given flavor 
susceptibility to the total one. Our analysis provides model-independent predictions for the mass, temperature and 
isospin dependence of those quantities, which should be recovered in lattice analysis as they approach the continuum 
limit. In accordance with the behaviour of the corresponding order parameters, the disconnected susceptibility shows 
a linear growth at low and moderate temperatures, infrared divergent near the chiral limit as T/M^r, whereas the 
connected one is infrared regular but survives the chiral limit as a growing T^/M^ behaviour. The chiral or infrared 
limit gives qualitatively the behaviour as the temperature approaches chiral restoration but numerically is not a good 
approximation for physical pion masses. The higher order isospin breaking corrections are quadratic in m„ — and 
are enhanced in the chiral limit for the disconnected susceptibility, as long as rus and niu/md remain flxed. The ChPT 
susceptibilities reproduce the growing T-dependence at low and moderate temperatures in a model-independent way. 
Although they do not show the peaks expected near the transition, our small mass analysis allows to infer that the 
disconnected part should have a more pronounced peak than the connected one, the latter expected to present a 
rather soft behaviour. This difference can be interpreted from the different order parameters that fluctuate in each 
case: the chiral quark condensate for the disconnected piece and the isospin-breaking one in the connected case. In 
the formal SU{2) limit nis ^ oo the connected contribution becomes temperature independent, like {uu — ddtjx- Our 
analysis for the susceptibilities is consistent with previous related work in the literature. 
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